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ABSTRACT 

•  * 

Conventional  (or  restricted)  Hartree-Fock  wave  functions  have 
been  obtained  for  Ga*.  Ge'*^''’,  Br*,  Rb+  and  the  neutral  atoms  Zn ,  Ga ,  Ge , 
As,  Se ,  Br  and  Kr.  Analytic  Hartree-Fock  methods  utilizing  a  version 
of  Nesbet's  symmetry  and  equivalence  restrictions  ,«were  used  to  obtain 
the  solutions  in  analytic  form  (suras  of  powers  of  r  times  exponentials),’ 
Results  are  given  in  the  form  of  wave  functions,  one-electron  energies, 
and  total  energies.  Comparison  is  made  with  the  earlier  results  for 
Zn  and  Kr,  the  only  atoms  in  this  group  for  which  H-F  solutions  exist. 
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•  *  * 

I .  INTRODUCTION 

« 

It  has  become  Increasingly  apparent  that  Hartree-Fock  wave  functions 
for  free  atoms  and  Ions  can  play  a  vital  role  in  providing  approximate 

m 

descriptions  of  molecular  and  solid  state  phenomenon.  A  striking  example 
of  this,  and  one  which  was  quite  unexpected,  is  the  successful  use***  of 
free  ion  Hartree-Fock  wave  functions  to  describe  the  recent  data  on  isomeric  * 
shifts  (l.e.,  total  s  electron  density  at  the  nucleus)  which  have  been 
determined  by  Mossbauer  experiments.  It  goes  without  saying  that  basic'to 
any  such  applications  is  not  only  the  existence  of  such  functions  but  their 
availability  in  a  form  suitable  for  computation.  In  some  of  our  own  studies 
involving  the  4p  atoms  the  complete  lack  of  H-F  wave  functions,  for  this  •> 
part  of  the  periodic  table,  was  only  too  obvious.  In  order  to  overcome  this 
difficulty,  we  have  determined  Hartree-Fock  solutions  for  Ga'*’,  Ge*"*^,  Br' ,  Rb"'' 

a 

and  the  neutral  atoms  Zn,  Ga ,  Ge,  As,  Se ,  Br  and  Kr.  of  these,  Hartree-Fcfcck 
(H-F)  results  have  been  previously  available  for  only"Zn'*’  and  Kr*®*. 

These  are  conventional  or  restricted  Hartree-Fock  solutions  In  that  one- 
electron  functions  of  the  same  shell  are  constrained  to  have  the  same  radial 
dependence. *^ *  - 

As  we  have  said,  aside  from  their  own  Inherent  interest  as  a  descrip¬ 
tion  of  the  electronic  structure  of  free  atoms,  a  major  purpose’  of  such 
calculations  is  to  supply  a  starting  point  for  further  investigations.  The 
spherical  ion  results  have  been  used  to  obtain  the  Sternheimer  quadrupole 
antishielding  factors*®*  (7oo’s)  as  these  extend  our  knowledge  of  theoretical 
values  of  to  ions  for  which  only  estimates  could  be  made  previously*** 
fexaapt  for  Rb^  for  which  Hartree  calculations  exist.)  The  H-F  results  have 
been  utillaed  in  saveral  other  investigations.  Atomic  scattering  factors  • 
have  been  obtalnad*'*^  and  the  Ge  results  have  supplied  a  starting  point  in 

O  • 

an  effort  to  improve  on  core  and  valence  electron  self-consistency  in  ortho- 
gonalized  plane  wave  calculations***  for  germanium. 
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Analytic  H-F  methods,*®’  utilizing  a  version  of  Nesbet’s  symmetry  and 
« 

equivalence  restrictions**’®’  have  been  used  to  obtain  the  wave  functions. 

9> 

»  » 

Some  details  of  the  method  are  disftussed  in  Section  II^  A  fullac  discus¬ 
sion  has  been  given  prev loosl^ ”  wher&  H-F  results  for  gp  shell  atoms  were 
reported.  The* H-F  results  appear  in  Section  III  along  with  comparisons  o£ 
the  prev  1 0U3 1 5*'  aval  labl  c  Zn  and  Kr  functions. 


II.  DESCRIPTION  OF  THE  CALCULATION 

Six  of  the  eleven  stOBs  and  ions,  for  which  H-F  results  will  be  re¬ 
ported,  are  closed  shell  lone.  These  can  be  straightforwardly  handled”®’ 
but  the  remaining  five  Ions  with  their  unfilled  4p  and  filled  2p  and  3p 
shells  raise  a  problea  of  maintaining  orthogonality  between  p  orbitals. 

We  have  discussed  this  natter  at  length  when  reporting’”’  3p  atom  solutions 
and  the  reader  Is  referred  there  for  a  discussion  of  this.  We  have  utilized 
the  sane  version  of  Nesbet's  symmetry  and  equivalence  restrictions*’®’  as 
was  used  In  the  3p  atom  investigation.”’’  This  procedure  does  not  yield 
the  Pest  possible  H-F  results  but  it  has  been  our  experience,”®’  and  that”^’ 
of  Investigators  who  have  used  the  numerical  H-F  method,  that  the  H-F  eigen¬ 
functions  are  Insensitive  to  the  way  in  which  the  orthogonality  "problem"  Is 
treated.  Because  of  this,  application  of  the  symmetry  and  equivalence  re¬ 
strictions  yields  results  which  are  only  negligibly  Inferior  to  the  "best 
results."  The  discrepancies  introduced  are  far  smaller  than  the  difference 
between  the  Hartree-Fock  and  the  "true"  nany-electron  wave  functions.  It 
is  our  opinion  that  if  one  requires  superior  many-eiectron  eigenfunctions  or 
eigenvalues  one  should  be  prepared  to  obtain  these  by  going  beyond  thecon- 
ventlcnai  Hartree-Fock  formalism. 

As  already  noted,  we  have  utilized  the  analytic  Hartree-Fock  method  to 
obtain  cuf  jgsulta.  This  method  msws  standard  matrlji  techniques  to  obtain 


orthonoriaal  anajyti®  Hafferce-F^offk  radial  «»rt)lt'ais,  of  th®  forrag 
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and  Cue  basis  functions  Rj's,  are  of  ttie  form: 


Rj  (r)  -  ^  e  (3) 

0 

wbere  t  Is  the  dne«electron  angular  momentua  quantum  number  appropriate  for 
the  one-electron  orbital  of  wliicb  Uj(r)  is  the  radial  part.  The  Nj  Is  a 
normalization  constant  and  is  expressible  in  terms  of  the  other  parameters, 
l.e. : 
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Uj(r)*s  of  comnon  4  value  are  constructed  from  a  common  set  of  nj(r)*s. 
Given  the  basis  sets,  l.e.  the  Rj(r)*s,  the  problem  is  reduced  to  solving 
the  Ilartree-Fock  integro-dlfferentlal  equations  for  the  eigenvectors  (the 
Cj^^’s)  and  their  eigenvalues,  liils  is  done  by  straightforward  matrix 
(lladonnl  izatlon  and  manipulation  and  avoids  the  problems  of  numerical 
accuracy  inherent  In  the  integrations  of  the  numerical  H-F  method. 


Recently ’  several  workers  have  overcome  the  problem  of  numer¬ 
ical  integration  accuracy.  In  particular  Wayers’^®’  has  not  only  obtained 
highly  accurate  numerical  11-F  solutions  but  he  has  also  accurately  evaluated 
(a  not  easy  task!)  the  one  and  two-electron  integrals  making  up  the  total 
energy.  This  will  allow  detailed  comparisons  tn  he  made  for  the  first  time 
between  numerical  and  analytic  results. 
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IQ,  general  the  twe  approeffhe^  aje  CfJmplementary ,  sometimes  on%  and 
^omstiresa  the  other  provlflas  a)  the  isore  ciiitable  method  of  computation 
snd/of  b)  resnlts  in  a  siore  convealent  form,  Vi'e  have  found  the  analytic 
approach  to  be  very  Hsefwl*^*’  wfien  going  beyond  the  conventional  free  atom 
H-P  fomallSB.  We  have  also  found  it  convenient  to  have  the  resultant  wave, 
functions  In  analytic  form. 

In  using  the  analytic  approach  we  have  replaced  the  problem  oT  accuracy 
of  nuQcilcal  Integration  by  the  problem  of  choosing  adequate  basis  sets. 

Let  vs  consider  this  loatter  now.  First  there  Is  the  question 
of  the  size  of  the  set.  A  seall  set  Is  desirable  because  of  economy  In 
computer  tine  and  retains  the  advantages  of  wave  functions  of  analytic  forn. 
These  advantages  cone  fron  the  ease,  accuracy,  and  convenience  with  which 
natrlx  elements  can  be  obtained  if  the  functions  are  In  analytic  form. 

Large  basis  sets  allow  greater  accuracy  of  solution  (providing  that  we  main¬ 
tain  sufficient  linear  Independence  among  the  basis  set,  otherwise  errors 
accumulate  during  matrix  dlagonalizatlon) .  Having  made  the  choice  of  the 
size  of  the  basis  set  there  is  then  the  problem  of  choosing  the  Individual 
Rj’s.  in  the  present  work  we  have  relied  heavily  on  earlier *^*  ^’’*  II-F  in¬ 
vestigations  for  the  Iron  series  ions  and  so  we  will  review  these  briefly. 
The  first  investigation'***  involved  obtaining  H-F  solutions  for  Iron  series 
Ions  In  Sd"  (l.e.  no  4s  electrons)  configurations.  The  basis  sets  were  the 
largest  that  could  be  fitted  onto  the  computer  that  was  used.*'®*  Their 
size  fixed,  they  were  obtained  by  an  extensive  series  of  computations  where 
the  basis  functions  were  varied.  This  search  concentrated  on  obtaining 
good  basis  functions  for  V  through  Nl.  This  investigation  was  followed  by 
a  series  of  II-F  calculations  for  the  neutral  4s^  Iron  series  atoms. 

These  were  done  on  a  cooputer  of  larger  capacity  and  the  basis  sets 
obtained  by  adding  additional  s-llke  Hj>s  (for  the  consbrVQUon  of  tb%  ouspr 
loop  of  the  43  orbitals)  to  the  othetbtise  nnmacllieiJ  bests  set^l  ei  the  ffrst 
Investigation.  An  extensive  bas^  s  function  search*was  thus  avoldeij  at  the^ 
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^■.Oa*  (ft  fn  tli<»  ,«li<}rttominR»  of  tiie  fif5t  efroffc.  Sntiireff!»eDt.l,7  m^re 

ECCurata  analytic  H-K  calculation#  w<j,re  (dene  fer  Zn  (3d^®4s^,  ^S)  whii^h  is  ^ 
reported  here  and  for  Fe  (3fl^4s2.5D) ,  Mn++  (sd®,  ®S)  a.n<3*Cu+  *S)  which appee* 

elsewhere.****  Th<»  greatest  dlscrepaacies  betweea  calculations  for  the  same 
state  occur  lor  Cu*  and  Zn  which  lie  outside  the  group  emphasized  In  the 
first  basis  function  search.  '■ 

The  present  computations  were  done  on  the  I.D.M.  104  computet  at  Avco 
whose  capacity  (32,000  words)  allowed  the  use  of  larger  basis  sets  and  we 
have  taken  advantage  of  this  capability.  In  fact,  the  choice  of 
Rj's  Is  based  on  the  Iron  series  Rj’s  and  on  an  extensive  basis  function 
search  for  Kr  followed  by  tests  on  the  resultant  sets  for  other  atoms.  As 
a  result  the  ll-F  functions  appearing  in  thl.s  paper  are  more  accurate  than 
the  earlier Iron  series  solutions.  The  parameters  defining  the  Rj*s 
appear  In  Table  I.  We  have  reported  (and  used)  screening  constants  (Zj's) 
with  four  «lglt.s  after  the  decimal  point.  This  does  not  nean  that  Zj's  were 
uniquely  established  to  this  many  digits.  The  Investigations  varying  the 
Zj's  carried  this  many  digits  and  since  these  were  kept  In  the  final  calcu¬ 
lations  they  are  reported  here. 

III.  RESULTS 

The  eigenvectors  (Cjj’s)  which  define  the  Uj(R)’s  In  terras  of  the 
Rj<r)’6  appear  in  Table  II.  Note  that  the  Cj^j’s  are  given  for  norma llzed 
Rj(r)*s.  The  C|j’s  have  not  been  uniquely  established  to  the  number  of 
digits  quoted  but  with  these  digits  they  provide  well-orthonorwa llzed,  well- 
defined  Hartree-Fock  orbitals.  Total  energies,  one-electron  energies  (Cj^'s), 
ore-electron  kinetic  plus  nuclear  potential  energies  (K{’s) ,  and  two-elcctron 
F*‘(4p,4p)  integrals*®*’  are  listed  in  Table  III.  The  quantities  appearing 
there  are  accurately  evaluated  for  the  functions  defined  In  Table  II.  TSie 
Nartree-Focls  one-glectron  and  total  energies  have  not  been  rigorously  ^ 
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•  established  to  the  number  of  digj^ts  quoted,  ^:r'^e  accurate  solution  of  the 

•  Hartree-Fock  equations  (e.g.  through  the  use  of  Hifgar  improved  basis  setsj 

will  affect  t^e  la^t  one  or  two-digits  of  th§5a  Ctiian^ities ,  iTote  that 

Hartrea  atorafa  unite  have  been  us6d»(l  ati  =  S  Ry) .  *  *  * 

»  * 

*■  ■% 

*Fhe  Kj’s  which  do  n®t  include  in,tet.clectronic  terras,  offer  a  Sensitive 

•  test  of  wave  function  behavior.  For  fixed  Z,  a  contraction  in  a  ,r)  causes 
an  increase  in  the  aagnltode  of  its  Kj^.  inspection  of  the  Kj*s  of  Table  III 
shews  in  general  that  the  U^*s  bf  any  one  element  contract  when  going  tc  a 
state  of  higher  positive  ionisation.  This  is,  of  course,  not  surprising. 

Of  greater  Interest  is  the  exception  to  this  rule  which  occurs  for  the  2p 
and  3p  orbitals  on  going  from  a  state  with  one  or  more  4p  orbitals  (such  as 
Ga  or  Qe)  to  a  state  with  none,  (such  as  Ga^  or  Qe^*^).  The  requirements  of 
obtaining  three  simultaneous  (hence  orthogonal)  p-llke  eigenfunctions  for 
the  former  state  and  the  partial  relaxation  of  this  for  the  latter  has  caused 
a  reverse  trend  in  K2p  and  K3p.  .  It  should  be  noted  that  the  resultant  2p 
and  3p  orbital  variations  are  small. 

Let  us  now  compare  the  present  results  with  those  previously  available 
for  Kr  and  Zn.  Worsley’s  numerically  obtained*®’  e^'s  for  Kr  appear  In 
Table  IV.  Kj's  and  the  total  energy  were  not  obtained  by  her.  We  see  that 
there  Is  substantial  agreement  between  the  two  sets  of  results.  More  detailed 
comparisons  can  be  made  with  the  earlier  obtained*®’  analytic  Zn  results 
where  Kj's,  e^’s,  and  the  total  energy  are  available.  These  appear  In  Table  V. 
Comparison  of  e^'s  show  larger  discrepancies  than  were  seen  for  Kr.  This  Is 
due  to  the  previously  discussed  deficiencies  in  the  basis  set  of  the  earlier 
calculation.  If  one  chose  to  inspect  the  to  learn  of  orbital  variation 

as  has  often  been  done  (out  of  necessity)  in  the  past,  one  would  conclude 
that  the  4s  orbitals  have  changed  least  (the  ^45*8  being  in  best  agree'* 
ment  between  calculations)  and  that  the  3d  orbital  of  the  present  calcula¬ 
tion  Is  the  more  contracted  (since  Its  £3^  is  more  negative).  On  the  other 
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hand  inspection  of  the  Kj's  shows  that  neither  conclusion  is  correct.  The 
4s  orbitals  show  the  greatest,  not  the  least,  modification  and  the  present 

m  m  • 

calculation's  3d  orbital  is  expanded  and  not  contracted  relative  to  the 
earlier  results.  In  addition  to  inspecting  IndlvldUEel  ’ s,  sums  '  '  of 

(which  also  appear  in  Table  V)  have  in  the  past  been  used  4:0  estimate  the 
variation  in  total  enecg^.  Suc^h  sums,suggest  that  there  has  been  a  1  au 
change  In  total  energy  when  in  reality  a  change  of  one  fiftieth  that  size 
has  occurred,  la  other  words  the  cae>slectroa  energies  supply  an  often 
misleading  "yamstlck"  to  wave  function  behavior. 

*  "  « 
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and  p.  4i;  No.  18,  October,  1956,  p.  4,  Solid-State  and  Molecular  Theory 
Group,  Massachusetts  institute  of  Technology,  Cambridge,  Massachusetts 
(unpublished)].  Recently  C.  C.  J.  Roothaan  [Revs.  Modern  Phys,  32,  179 
(1960)1  has  extended  his  formalism  to  cover  the  nonclosed  shell  case  for 
th'e  conventional  restricted  Hartree-Foclc  method  where  nonzero  off-dlagonai 
Lagrange  multipliers  occur.  S.  Huzlnaga  (Phys,  Rev.  120,  866  (I960); 
ibid  122 ,  131  (1961))  has  extended  Boothaan's  formallsB  further. 

^®R.  K.  Nesbet,  Proc.  Roy.  Soc.  A230,  312  (1955). 

E.  Watson  and  A.  J,  Freeman,  Phys.  Hev.  (to  be  published);  Section  II 
and  the  appendix  of  this  paper  supply  sufficient  Information  for  the  reader 
to  construct  the  Hartree-Fock  equations  solved  in  the  present  paper. 

’^See  D.  R.  Hartree,  The  Calculation  of  Atomic  Structures,  (John  Wiley  a 
Sons  Inc.,  New  York,  1957)  for  details  concerning  the  conventional 
Hartree-Fock  formalism  and  for  the  derivation  of  Hartree-Fock  equations.® 


^^Thls  observation,  is  based'on  fche«wj>fTt  reported  !!n  Pef.  11  and  In 
R.  k:  Nesbet'and  R.  E.  Watson,  Ann.  Phys.  (New’  YorSJ  2^0  j;i960). 

• 

^“^In  two  cases  (JD!  R^.  Hartreft  and  W.  Hartfes,  Pfaft.  Roy.  Soc.  A193 .  599 
(1948)*,  W.  Hartree*,  D.  R.  Hartree  and  M,  Manning,  pfiys.  Rev.  857  »> 
(1948),  small  "off-diagonal"  .Lagrange  multipliers  were  Included  for  " 

orthogonality,  otherwise  [d.  R.  Hartree  and  W,  Hartree,  Proc,  Roy.  Soc. 

A164 ,  167  (1938,);  D.  A.  Goodings  (to  be  published);  and  unpublished  work 
of  D.  F.  Mayers*  (1958)  the  problem  was 'Ignored ,  le.,  the  multipliers 
were  set  equal  to  zero  and  orthogonality  was  reasonably  maintained. 

F.  Mayers  (to  be  published)  has  obtained,  among  others,  numerical 
Hartree-Foc!?  results  for  Zn  but  at  the  time  ol  writing  these  results  ,, 
were  in  a  preliminary  stage  and  thus  were  not  available  £or  comparison, 

1®D.  A,  Goodings  (to  be  published)  has  obtained  accurate  numerical 
"unrestricted"  Hartree-Fock  solutions  for  a  number  of  low  2  atoms. 

^"^For  example,  the  analytic  method  was  used  to  obtain  Hartree-Pock  solu¬ 
tions  (1)  for  ions  in  external  "crystalline*  Helds  [R.  B.  Watson,  Phys. 
Rev.  U7 ,  742  (1960)  and  R.  E.  Watson  and  A.  J.  Freeman  Phys.  Rev.  120 , 

1134  (1960).];  (2)  for  a  non-magnetlc  ion  In  the  "exchange  field"  of 
neighboring  magnetic  ions  [A.  j.  Freeman  and  R.  B,  Watson,  Bull,  Am.  Phys. 
Soc.  ^  234  (1931)];  and  (3)  to  obtain  one-electron  orbitals  which  are 
not  separable  into  a  product  of  a  radial  and  an  angular  function 
C.  Sonnenscheln ,  (to  be  published),  in  (2)  and  (3)  the  analytic  method 
has  several  advantages  In  making  the  computations  possible. 

1®R.  B.  Watson,  Phys.  Rev.  US,  1634  (1960). 

Whirlwind  I  was  used. 

2°For  definitions  see  B.  U.  Condon  and  G.  H.  Shortley,  The  Theory  of  Atomic 
Spectra,  (Caabrldge  University  Press,  Cambridge,  1953),  p.  177. 

^^The  Mn*"*  results  appear  in  an  investigation  [R.  B.  Watson  and  A.  J.  Freeman 
(to  be  published)]  of  Uyperflne  effects,  those  for  Cu*  In  an  investigation 
of  quadrupole  polarizabilities  [A.  j.  Freeman  and  R.  B.  Watson,  Bull.  Am. 
Phys.  Soc.  ^  168  (1931)],  and  the  Fe  results  ate  unpublished. 
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Table  IV  Kr  one-electron  energies  for  the  present  calculation 
and  as  obtained  by  Worsley^  «,  ^ 


0 

Worsley's  results®  #, 

present  calculation 

0 

^IS 

-520*  au 

-520.167  au 

• 

^2S 

^  -  69.9  ,  " 

-  69.904 

0  • 

*3S 

-  10.85 

-  10.850 

• 

i 

f  0 

®48 

“  1  •  1 51 5 

-  1.153 

^2p  . 

-  63.1 

-  63.011 

^3P 

-  8.335  * 

-  8.332 

«4p 

-  0.53 

'  -  0.524 

CO 

-  3.80 

-  3.826 

®See  ref.  3 
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Table  V  Zi»  one-olectroa  energies  one-electron  +  nuclear 

*  potential  energy  Integrals  {ifi'sj,  sums  of  ofle-elu^tron’energles, 
and  total  energies  for  the  present  and  earl  ler^eagali^tlc  Hartree- 
Fock  calculations.  * 


SarSler  cateuiation^ 


Present  cal^ifl  atTofi 


-353. 2S1.  a« 


-353. 299  au 


-  in. 319 


-  lilt.  358 


-  5. 600 


-  5.635 


-  0.266 


-  0.  29  f 


38.882 


38.921 


-  3.804 


-  3. 837 


-449. 849 


-449.848 


-1 10.  569 


•110.571 


-  43.  136 


43.  155 


12. 060 


-  12.218 


-109.912 


-108.913 


-  40.996 


-  40.881 


-  34.842 


-  34.722 


Total  Energy 


-1777.823 


-1777.843 


'1070. 568 


-1071.514 


®8ee  ref.  2 
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